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1 Introduction 

1 derive the third order optics of quadrupoles. I transform away the deriva¬ 
tives of the strength function, and thus demonstrate that third order aber¬ 
rations are insensitive to fringe field shaping. The results can be used for 
efficient tracking to third order or for simple and quick evaluation of nonlin¬ 
ear effects. 

This is an extension to my work of 1997[[1], which covered both electro¬ 
static and magnetic quads, but only for the non-relativistic case. 

2 Theory 

A particle of charge q and mass m has canonical pairs of coordinates ((x, P x ), (y , P y ), (z , P z ), (t, E )). 
The equation coupling these is 

(.E — g$) 2 — | P — qA\ 2 c 2 = m 2 c 4 . (1) 

Here, $ and A are respectively the scalar and vector potentials; both are 
functions of ( x,y,z ), and both are zero on the reference particle’s orbit. 

Let us use units for time, energy and potential that sets resp. c = 1, 
m — 1, q — 1. This is permissible as long as we consider no processes that 
change mass or charge. 

(P-$) 2 -|P-A| 2 = l. (2) 

To find the Hamiltonian El = E, solve for E: 

H(x,P x ,y,P y ,z,P z -t) = <L+ \Jl + \P - A\ 2 (3) 
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2.1 In Frenet-Serret Frame 


The reference orbit is assumed to be in one plane and is generally curved 
with curvature h = 1/p. The transformation is conventionally made to the 
Frenet-Serret coordinate system where the longitudinal coordinate s is in 
the reference orbit direction, so h = h(s), x is radially outward, and y is 
perpendicular to the bend plane. Ruth [2] shows that then the Hamiltonian 
is 

H(x, P x , y, P y , s, P s] t) = $ + \jl + {P x - A x y + {P y - Ayf + ~ ^ )2 

(4) 

As conventional in beam and accelerator physics, we use the longitudinal 
coordinate s as independent variable. Then the Hamiltonian is — P s : 

H(x,P x ,y,P y ,t,E-s) = -A a -(l+hx)y/(E - <f>) 2 - 1 - (P x - A x ) 2 - {P y - A y ) 2 

(5) 


2.2 In Differential Coordinates 


The “reference particle” has x — y — 0 and P x = P y = 0. In the following, we 
use the traditional symbols /3, 7 and hence also /3y for the reference particle’s 
speed, energy, and momentum, respectively. 

The Hamiltonian is awkward because it mixes small dynamic quantities 
x,y,P x ,P y with a large one E. We only care about particles with a small 
A E deviation from the reference energy 7 , and a small At deviation from the 
reference time to = S /P- We do this with a canonical transformation from 
(t, —E) to (At, — A E). The generating function is 


F(t, —AE) — t-- (-A£- 7 ) 


The new Hamiltonian is 

H, = IIs 


dF 

ds 


(3 


= H„ 


7 + AE 


( 6 ) 

(7) 


Furthermore, we introduce new coordinates (r, P T ) in place of (At, AE), with 
r = /3At, P T = AE/fi. This results in a new “time” coordinate r being the 
distance ahead of the reference particle, and the “energy” coordinate being 
the momentum deviation w.r.t. the reference particle. 
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The Hamiltonian is then 

H s (x, P x ,y, P y , t,P t ;s) = 


( 8 ) 


P t -A s -{ 1 + hx)yj{ 7 + (3P T - $)2 - 1 - Jp x - A x ) 2 - (Py - A y y 


2.3 Re-Normalize Momenta 


We now change the units of momentum to [3 7 . This will have the advantage 
that outside of the regions of electric and magnetic holds, we have P x = x ', 
P y = y', P r = AP/P where primes are derivatives w.r.t. s. Further, we 
rescale scalar potential by a factor /3 2 y, and vector potential by a factor ^ 7 . 
The result is 


H = P T -A s -(l+hx)J Q + P(P T - $)) - ~ - (P x - A x y 
or 

H = P T -A s -(l+hx) yjl + 2 (P T - $) + /3 2 (P r - <t >) 2 - (P x - A x ) 
This is the general Hamiltonian. It is exact. 


- (Py ~ AA 

(9) 


2 - (Py - A ,) 2 
( 10 ) 


2.4 Potentials’ Scales 


For handy reference, here are the definitions of the scaled potentials in terms 
of the unsealed (subscript u): 


$ = 
A = 


q 


/ 3 2 'ymc 2 


$u (x,y,s), 


Q 


(3ymc 


A u (x,y,s) 
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2.5 Relativistic Limits 


This H also has the nice feature that the non-relativistic and ultra-relativistic 
limits are simple: 

P « 1 : H = P T - - (1 + hx)^J 1 + 2 (P T - $) - (P, - A x ) 2 - {P y - A y ) 2 

7 > 1 : H = P T - - (1 + hx)yj{ 1 + P T - $) 2 - (P x - A x ) 2 - (P y - A y ) 2 

( 12 ) 
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2.6 Example: Field-free, curvature-free 


H s = P r - ^1 + 2 P r+ ^p2_p2_p2 
Then P' x = P' y = P' T = 0, and to first order: 

^ = dH_ = _ ^ _ } 

dP x ^/1 + 2 P T + p*P? -Pl-P>~ 

and similar for y', 

^, = 1 _ 1 +/3 2 Pr _ _ K 

^1 + 2 P T + /3 2 P? - PI - PI ~ 7 2 


( 13 ) 


(14) 


(15) 


2.7 Straight elements, ignore longitudinal 

Let us now confine ourselves to straight elements ( h = 0) and concentrate 
only on transverse. Then for magnetic elements, we have 

H — —A s — yjl - {P x - A x y - ( P y - Ay) 2 (16) 

and for electrostatic elements, we have 

H = - S J 1 - 2$ + /W - PI - PI (17) 


3 Electrostatic Quads 

Compared with equation 2 of the 1997 paper|T], we notice an extra term 
/3 2 $ 2 in the square root of eqn. cfl 

We expand the square root to 4th order in coordinates and ignore the 
constant: 

H « ^(2$ - /3 2 $ 2 + Pi + P 2 y ) + i(2«h + Pi + Pi) 2 . (18) 

To the same order, Laplace’s equation gives for the expansion of the 
quadrupole potential: 

* = as) 

Where is a also a factor of 2 because the scaling of eqn. [lT] differs from the 1997 scaling 
by this factor 
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The expanded Hamiltonian, correct to 4 th order is 


H 


k(x 2 - y 2 ) P 3 


P 2 


+ -f + -f+ 


, (P* + Py 2 ) 2 , M ^- l / 2 )(^ 2 + n 2 ) 

+ -^ + -4-“ 

k 2 (x 2 -y 2 ) 2 k"(x A -y A ) 


8y 2 


24 


( 20 ) 


3.1 k" is a “Fringe Field Effect”? 

The trouble with applying this to simple cases like thin lenses and hard-edge 
limits is the presence of k"(s), which becomes singular in those limits. In 
most cases, one sacrifices physical insight and simply traces particles with 
this Hamiltonian, using a more-or-less realistic function k(s). For example, 
the approach taken in GI0S[3| is to leave it up to the user to specify ‘fringe 
field integrals’ such as f k 2 ds through the fringe fields. However, this leaves 
one quite vulnerable to error; different integrals may not be realistic or con¬ 
sistent with each other. Moreover, if one needs to solve Laplace’s equation 
to find fringe field integrals, one might as well use the solution directly in 
a ray-tracing code. If one does go through this exercise, one discovers that 
the higher order aberrations are relatively insensitive to the ‘hardness’ of the 
quadrupole edges. This leads one to suspect that the aberrations are domi¬ 
nated by an intrinsic effect which has nothing to do with the detailed shape 
of the fringing field. Such is indeed the case. 


3.2 k" can be transformed out! 

It turns out to be possible to find a canonical transformation which eliminates 
the derivatives of k(s). In our case, we wish to retain terms to 4 th order in 
the Hamiltonian (3 rd order on force), and the transformation (x, P x , y, P y ) —>■ 
(X, Px, Y, Py) has generating function 

G(x, Px, y, P Y ) = xPx + yPy + - y 4 ) - hx 3 P x - y 3 Py). (21) 

24 o 
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To the same order, this yields the transformation 

x = X + \x 3 

k k' < 22 > 

P* = Py - -A 2 P* + -X 3 . 

2 o 

The //-transformation is obtained by replacing x, P x , A", P x with y , P yi Y, P Y 
and k with —k. Note that outside the quadrupole, the transformed coordi¬ 
nates are the same as the original ones. 

This yields the transformed Hamiltonian H*: 


H* 


= |(A' 2 - Y 3 ) + l -(P 3 + P 3 )+ 

+ \‘X 2 x + P'i-f-- i (X 1 + Y 2 )(Pl- 

(7 —3/? 2 )fc 2 4 4 , (1 — P 2 )k 

+-24- (A +i } - r~ 


Py)+ 

-X 2 Y 2 . 


(23) 


We can identify the terms: the first two are the usual linear ones; the third 
term is not related to the electric held (it is small and due to the fact that 
x'P x or, equivalently, tan# ^ sin#); the 4 th term is also small and arises 
because a particle going through the quadrupole at an angle is inside the 
quad for slightly longer than one which remains on axis. See ref. (4] for more 
complete physical derivation of the individual terms. 


3.3 Thin lens, Hard Edge Formulae 


The dominating higher order terms are the last two terms in eqn.[23l Since 
there are no derivatives of k , we can directly write down the aberrations in 
the thin-lens limit: 


A P x 


-1 {7 — 3 /? 2 3 
PL { 6 X 



(24) 


with a similar expression for A P y . L and / are the quadrupole’s length and 
focal length. (Actually, it is more accurate to replace jq- with f k 2 ds.) 

The fractional focal error is found by dividing by the linear part Ao P x = 
-x/f: 


A /, _ 1 / 7 — 3J3 2 2 
/ fL\ 6 

for x, and similarly for y. 



(25) 
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3.4 Physical Interpretation: Speed Effect 

It is interesting and instructive to deconstruct the final result [24] to derive 
physical origins for these terms. We do this in the thin lens limit. 

Referring to the untransformed Hamiltonian (l20|h we can identify the 
term with r y~ 2 as due to a “velocity-gain” effect: particles entering the electric 
field have their speed changed because of longitudinal field, for example, those 
entering near the like-charged electrode are slowed so spend a longer than 
normal time in the focusing field. This effect disappears in ultra-relativistic 
limit. The contribution to P x from this effect in thin lens limit is 

aau = -hTxfi! j (26) 

leaving only |x 3 inside the parentheses of eqn.ET to account for. 


3.5 Physical Interpretation: k" Effect 

The direct effect of the k" term in the potential and the Hamiltonian (1201) 
can be found from integrating by parts: 


AP x \k" — 


1 

2 


dH 

dx 


ds 

k" 


kx 2 x"ds ~ 


1 

6 

1 

2 


k"x 3 ds = 


k 2 x 3 ds ~ 



(27) 


3.6 Physical Interpretation: Ax Effect 

The remainder is now — ^ f k 2 ds. This originates from a small, subtle and 
largely overlooked effect; I overlooked it in my earlier workjTj. It originates 
from a shift in x experienced in the quad. 

From [27] above, we also find 


x 


n 


k"x 3 

~ 6 ~ 


which in thin lens approx can be integrated directly to obtain 


Ax = 


kx 3 

IT 


(28) 


(29) 
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Another way to see this is from eqn.|52J Since the transformed variable 
X does not see any shifts due to derivatives of k, it is unaffected on passing 
through the fringe held. But if A" is not shifted, then x must be shifted by 
kx 3 / 6. 

As x is shifted, there results a different overall focus effect A P' x = kAx: 



(30) 


Equations 1261 [2711301 when summed, givel24l Q.E.D. 

4 Magnetic Quads, Scalar Potential 

We can use the same scalar potential for magnetic as for electrostatic (flUjl . 
but rotated by 7t/4: 



(31) 


To find the vector potential, we follow Venturini-Abell-Dragt[5] (VAD) and 
express it first in polar coordinates: 

T (r, 6, s ) = Qr 2 - ^-r 4 ^ sin 2 0 (32) 

4.1 Vector Potential, Gauge Choice 

But instead of VAD’s gauge condition A e = 0, we set A r = 0. Then 


dA s 1 <9T 

dr r 86 

1 d(rAg ) 9T 


(33) 


r dr ds 


and we can find A simply by integrating. This results in 


Ag = —r 3 sin 26 
8 


k' 



(34) 
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or in Cartesian: 


y y u yt 

A I ™ 2 h 2 ™ / 2 2\ , h f A 43 

4 = I -yv . y y> -y 1 -v) + - v ) 


(35) 


It is interesting to compare this with the VAD[5j vector potential 
A = (T y - xy>) , J y-y - y 3 ) , y - » 2 ) + h y - a 4 )) (36) 

The two vector potentials differ by the gradient of the following function: 


\(x, y, s) = — y - y 4 ) k' 


(37) 


4.2 Lorenz Gauge 

Neither of these two vector potentials satisfy V • A = 0, the Lorenz gauge. 
We can add any multiple of x and it turns out that adding Vy/3 to eqn.l35l 
does satisfy Lorenz gauge to required order: 


A = 



(38) 


4.3 Magnetic Hamiltonian 

I choose to use eqn.l35l because it’s the simplest. Then the Hamiltonian can 
be written: 


H = 


k(x 2 ~y 2 ) , Pi , P y 


2 

(P^ + Py) 2 k'xy(yP x -xP y ) k"(x A -y A ) 
+ 8 + 4 48 

The generating function which will eliminate derivatives of k is 

G(x, P x , y, Py) =xP x + yPy + “ V*) + 


+ — + —+ 
2 2 


k_ 

12 


[(x 3 + 3 xy 2 )P x ~ (3 x 2 y + y 3 )P Y ] , 


(39) 


(40) 


9 






which, to the same order yields transformation 


x = X + — (X 3 + 3 AT 2 ) 

U k k> (41) 

P x = Px~~ [( X 2 + Y 2 )P X - 2 XYP y \ + 

and similarly for (y, P y ). The transformed Hamiltonian is 
h - =|(A' 2 - y 2 ) + 1(P 2 + p 2 )+ 

+ \(Pl + Pi? - |(a 2 + y 2 )(pf - Pi) (42) 

+ -(x 4 + Y 4 ) + ^x 2 y 2 . 


Notice the similarity to eqn.[23j in fact all terms are identical except the last 
two, which only differ in their coefficients. 

Applying the same procedure as in the electrostatic case, we find 


A P x = 




Or the fractional change in focusing strength: 


A/ x 

/ 


1 

Jl 



(43) 


(44) 


where L is the effective length. 
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